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Abstract 

The generation of harmonics by atoms or ions in a two-color, coplanar field configuration with 
commensurate frequencies is investigated through both, an analytical calculation based on the 
Lewenstein model and the numerical ab initio solution of the time-dependent Schrodinger equation 
of a two-dimensional model ion. Through the analytical model, selection rules for the harmonic 
orders in this field configuration, a generalized cut-off for the harmonic spectra, and an integral 
expression for the harmonic dipole strength is provided. The numerical results are employed to test 
the predictions of the analytical model. The scaling of the cut-off as a function of both, one of the 
laser intensities and frequency ratio n, as well as entire spectra for different w and laser intensities 
are presented and analyzed. The theoretical cut-off is found to be an upper limit for the numerical 
results. Other discrepancies between analytical model and numerical results are clarified by taking 
into account the probabilities of the absorption processes involved. 
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I. INTRODUCTION 



The possibility of obtaining high frequency radiation through the interaction of a laser 
field and an atom is a topic that has been extensively addressed during the last two decades 
from both a theoretical and an experimental point of view (see [|IJ for recent reviews). 
Harmonic generation from sources other than atoms, like linear molecules , ring molecules 
(e.g. benzene) || |J, nanotubes 0, and plasmas || have been also investigated. 

From the invariance of the Hamiltonian under dynamical symmetry operations, selection 
rules for harmonic generation can be elegantly derived @, 0- It turns out that apparently 
very different target and field configurations yield selection rules of the same type. Let us 
consider, e.g., a ring molecule with N ions (e.g., N = 6 in the case of benzene) in a circularly 
polarized laser field of frequency uj. The electric field vector lies in the plane that is spanned 
by the molecule and which we parameterize through the polar coordinates p and (p. The 
corresponding Hamiltonian (with the laser interaction taken in dipole approximation) is 
invariant under the dynamical symmetry operation 

A ( 2jt 2f \ 

P N =[p^P^^ V + -,t^—). (1) 

From this invariance follows 0, ^] that only harmonics of order 

n = gN±l, g eM+ (2) 

can be emitted. The harmonic radiation is circularly polarized and subsequent harmonics 
are alternately clockwise and counter-clockwise polarized. More complicated selection rules 
arise when also excited states are taken into account |7[]. Let us now turn to the actual 
target and field configuration examined in the present paper, namely the situation of an 
atom (or ion) in a circularly polarized two-color laser field of frequencies u and r]uj with 77 a 
positive integer number, and coplanar polarizations. In the case of counter-rotating electric 
field vectors the Hamiltonian is invariant under the same symmetry operation (P with TV 
replaced by rj + 1. The polar coordinates p and y? are with respect to the polarization plane 
now. In the case of co- rotating electric field vectors N has to be replaced by r] — 1. 

An appealing feature of the selection rule (|2|) is the fact that with increasing N, i.e., 
number of ions or frequency ratio of the two laser fields, respectively, less harmonics are 
emitted within a fixed frequency interval. This filtering effect may be accompanied with more 
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efficient emission of harmonics at short wavelengths which are of interest in spectroscopic 
applications, for instance. 

Harmonic generation in two-color fields has been studied both experimentally || and 
theoretically |g 0,[TT[. 

The present paper is organized as follows: In Sec. [TT] the theoretical modeling proposed 
by Lewenstein et al. ]T2|| for harmonic generation in the case of an atom interacting with a 
single linearly polarized field is extended to the two-color configuration, and several expected 
features of the harmonic spectra are deduced, among them the scaling of the cut-off and the 



dependence of the relative dipole strengths within a certain harmonic couple g. In Sec. |TJ ab 
initio numerical results obtained through the integration of the time-dependent Schrodinger 
equation for a two-dimensional model ion are presented and compared with the predictions 
by the analytical model. Finally, a conclusion is given in Sec. [1^. 
Atomic units (a.u.) are used throughout the paper. 

II. ANALYTICAL THEORY 

A theory of harmonic generation should answer mainly two fundamental questions: (i) 
which harmonics are emitted and (ii) which is the intensity of the harmonics as a function 
of the laser and target parameters. These questions have been addressed in [|12| for the 
case of an atom (in the single-electron approximation) interacting with a linearly polarized 
laser field (in dipole approximation). In a similar approach, the more general and more 
complicated case of elliptical polarization was studied in (single color). The elliptically 



polarized two-color field was addressed in |14| but the discussion of the cut-off law as well as 
the presentation of the numerical results were restricted to linear polarization there. Here, we 
focus on the case of two laser fields with circular polarizations and arbitrary integer frequency 
ratios and compare carefully the model predictions with ab initio numerical simulations. 

The electric field caused by the two lasers of frequency lo\ — u and u 2 = t\uj is assumed 
to be 

E{t) = (E x (t),E y (t),E z (t)) = (#! cos(wt) + E 2 cos(r]iut),E 1 sm(ut) - E 2 sm(rjujt),0) (3) 

where E\ and E 2 are the amplitude of the first and the second laser field, respectively, and 
the dipole approximation is applied. The two fields are oppositely polarized and coplanar. 
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The case of co-rotating electric field vectors will be discussed later-on. The vector potential 
A{t) = -J t E{t') dt' reads 

(EE E E \ 
— sin(atf) H sin( nut), cosiujt) H cosinojt), . 
UJ T]UJ UJ T]UJ J 

(4) 

Our starting point is the dipole moment along the direction n as it is calculated in the 



Lewenstein model (cf. Eq. (10) of Ref. [fP2| ) 

x H (t) = ij d£ J d 3 pn ■ d*(p - A(t))E(t') ■ d(p- A(t'))exp[-iS(p,t,t')] +c.c. (5) 

Here, 

dip) = (p\m, (6) 

and the action S(p, t, t') is given by 



In order to arrive at expression fl5|) several assumptions have been made in [12]: (i) among 
the bound states only the ground state plays a role in the evolution of the system; (ii) the 
depletion of the ground state can be neglected; (iii) in the continuum V(r) plays no role and 
the electron is treated like a free particle and can be therefore described py plane waves \p); 
and (iv) contributions from continuum-continuum transitions to harmonic generation can 
be neglected. 

For the field (||), (H), the general expression (||) evaluated for, e.g., n = e x , reads 

x(t) = i f dt 1 J d 3 p (E x {t')d x {p- A{t')) + E y (t')d y (p- A(t')) d* x (p- A{t)) exp[-iS(p, t, t')} 

(8) 

(the "+ c.c." is suppressed from now on). 

The integration over p is performed approximately by means of the saddle-point method, 
assuming that the major contribution to the integral is given by stationary points of the 
classical action, i.e, the points p^, p S y that satisfy 

%S(p 1 t,f) = 0. (9) 

One finds 

p x \t,T) = ^-( cos(cut) - cos(u(t - t))) + -^-( cos(7]ut) - cos(r]Uj(t - t))\ (10) 

UJ Z T \ J 7] Z UJ Z T \ / 



Py{t,r) = ^-(sm(ut) -sm(u(t-T))) - -^-f sin(^) - sm( V uj(t - r)j) (11) 

y LU Z T \ J 7] A UJ A r \ J 

where r = t — t' is the electron's "travel time." Introducing the stationary action 

S st (t, r) = S(p st , t, t - r) = ^ |M!bW + /p j df (12) 

with p st = (p S x,p S y, 0) we obtain, after the saddle-point integration over p, 

/•« / 7T \ 3/2 

x(t) = i / rfr — — - exp (-iS st (t, r ) r) - A t) 

jo \e + iT/2J 

x ^(t - r)d x (p s t(t, r) - A(t - r)) + - T)d y (p st (t, t) - A(t - r))) . (13) 



The factor with infinitesimal e in (p~3| ) comes from the regularized Gaussian integration over 
p. It expresses quantum diffusion of the released wave packet and damps away contributions 
from times r much larger than a laser cycle, allowing for the extension of the r-integration 
to infinity Pf. 



The stationary action ( 12 ) can be written in the form 



Setit, t) = C (r) + d(r) cos ((77 + l)u (t - r/2)) (14) 



where Cq(t) and C\{r) are given by 



^'-2^ (rV - 2 + 2c ° s( ™ )) + 2^r 



C (r) = 7 p r + — ^(r 2 ^ -2 + 2 cos(™)) + - 2 , (r? 2 rV -2 + 2 cos(??Ta;)), (15) 



Ck(r) = ^f 2 - sin (^r/2) sin(^r/2) - 2 ^f 2 3 sin ((77 + ]>r/2) . (16) 
r) z Ld^T 77(77 + l)ur 



The expression fll4|) for the quasi-classical action is very useful and interesting. The time 
dependence is given by just one term and through only one effective frequency which is 
(77 + l)u. This is consistent with the selection rule g(rj + 1) ± 1 obtained previously. Taking 
77 = 1 and E\ = E2 = E/2 the coefficients Cq(t) and Ci(r) calculated in \T2\ for the single 
linearly polarized field are easily recovered. 

As the semi-classical action is the integral over time of the kinetic energy, an expression 
for the energy gain of the electron is obtained by deriving S st (t, r) with respect to t, 

AE kin (t, r) = E kin (t) - EUt ~r) = dSst ^ r) 

= -(r; + 1)uC x {t) sin ( (77 + IV (^-)) ■ (17) 
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The value of the maximum energy gain is equal to the maximum of the function Ci(t) = 
(?7 + l)a;|Ci(r)|. Note, that Ci(r) depends on the product of the two electric field amplitudes 
only. If one of the two fields vanishes the resulting electric field is a pure circularly polarized 
field and, as expected, there is no possible energy gain and therefore no harmonics are 
emitted. 

We introduce 

U = ^ (18) 

which, for linear polarization and E\ = Ei — E/2, yields the well-known ponderomotive 
potential U p = E 2 /Au 2 . In Fig. [| the function C\(t)/U is plotted for 77 = 1,2,3,4. Writing 
the maximum energy gain as Ai?£^ x = j^U one obtains for 77 between 1 and 5 the values 
71 = 3.17, 72 = 1-28, 73 = 0.91, 7 4 = 0.67, 75 = 0.52. 

By expanding a part of the integrant in flT3]) in Fourier components, 



d*Mt(t, r) - A(t)) [E x (t - r)d x (p st (t, r) - A(t - r)) + E y (t - r)d y (p st (t, r) - A(t - r)) 

= ^6 M (r)exp(-iM^), (19) 



the dipole projection x(t) becomes 

r°° ( n \ 3/2 
x(t) = iE / dT \ —-7^ h ^ T ) exp(-iMa;t) 



xexp[-iC (r)]exp Ci(r) cos[(t] + l)u(t - r/2)]J . (20) 

In the case of an infinite laser pulse the coefficients 6m (t) are non-zero only for M = 
— 1 + A, 1 + A, —7] + A and rj + A, with A = m{r\ + 1) and m G M . This can be easily 
obtained writing down the dipole d(p) for a hydrogenlike ion ]nj. Moreover, one can see 
that |6m(t)| decreases very rapidly with increasing m and the leading terms are those with 
m — 0. Making use of a few variable changes and performing a Fourier transformation 
xk = i f% u / u dtx(t) exp(iKut) yields 



- ^S/;*( 7T ^) 3/2 ^-p( 1 (A--MW2)exp [ - i C„ W ] 



ir/2/ 77 

J w dtexp ^i ( — — - ^ ootj exp C\(t) cos(c<jt)^ . (21) 
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The integration over t in zero unless (K — M) /(rj + 1) is an integer number. In this latter 
case with the help of the Bessel functions of integer order J n (z) and taking into account that 



one obtains 



r e i,cos(0) e in0 rf = 2TTi n J n ( Z ) 

J —TT 



r°° ( 7T \ 3/2 

Zk = i£ / dT \ -77; b M {r) exp[i(K - M)ur/2] 

M Jo \ e + it/2 J 



x exv[-iC n (T)}i(—)jK-M [d (r)]. (22) 

17+1 

Given (X — M) / (77 -(- 1) — g with g £ J\f and considering the possible values of M one finds 
that xk is different from zero for K — g(r] + 1) + 1 and /<" = (7(7/ + 1) — 1. In the first case 
the contributing terms are those with M = 1 + A and M = —rj + A while in the second 
case they are those with M = — 1 + A and M = r\ + A. Hence, 

r f tv \ 3/2 

x g ( v+1)+1 = i 9+1 J o dr ( e + ir / 2 J exp[-iC (r)]exp[i^(r/ + l)^r/2] 

x ^(rJJ^dtr)] + i exp[i(r/ + l)^r/2]6_,(r) J fl+1 [d(r)]) , (23) 

f°° ( TT \ 3/2 

x 9( ^+i)_i = i 9+1 y o rfr I 6 + i r /2 ) ex P[- iC< o(7-)]exp[i^(r/ + l)wr/2] 

x (ft_i(r) J fl [Ci(r)] - i exp[i(r/ + l)wr/2]6,(r) J fl _i[Ci(r)]) . (24) 

By virtue of Eqs. (|23|) and ( p4|) it is seen that the selection rule fl2|) with ]V = 77 + 1 is 
automatically recovered. Note, that fl2"3"|) and fl2"4"|) are made up of two terms, one common 
(apart for coefficients which are the complex conjugated of each other) proportional to J g , 
and one proportional to J g+ \ for ( |2"3"D and to J 9 _i for (P4"|). The coefficients 6jw(r), to be 
calculated from ([Bp, are functions of the two laser fields E\ and Ei- In the limit E\ <C £2 
one has &±i(t) <C b± v (r) while ^> -E 2 implies fe±i(r) ^> ^(r). It follows that the 
intensities of the two harmonics g{r] + 1) ± 1 in a couple g are expected to be different for 
Ei <C E2 while they converge to the same value for E\ S> E2. 

III. NUMERICAL SIMULATIONS 

The numerical simulations were performed by integrating the time-dependent Schrodinger 
equation on a two-dimensional (2D) grid. Reducing the grid to 2D allows to run simula- 
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tions quickly on every modern PC and does not introduce qualitative modifications to the 
phenomena we are interested in here. Three-dimensional (3D) simulations are feasible but 
significantly more demanding. Related work on circular two-color stabilization of H in full 
3D has been published recently |TB ]. 



In polar coordinates (p,ip), length gauge, and dipole approximation the time- dependent 
Schrodinger equation under study reads 



~ TpTp ~ 2p^d^ ~~ d? + at ^ 

+ sm 2 (Qt)(Eipcos((f — tot) + E 2 pcos (</? + rjut) 



V(p,ip,t) (25) 



where the two laser pulses have a duration T = n/Q and a sine-square shape. Kt(p) is a 
"soft-core" 2D potential given by 

VM - (26) 

The parameters a and (3 can be tuned in order to adjust the ionization energy and the 
"smoothness" of the potential. In our simulations we used a = 2.44 and (3 = 0.20. These 
values provide an ionization potential of I p = 2.0, i.e., the one of real He + . The fundamental 
laser frequency uj was chosen 0.027T and the pulse length was T = 12600, corresponding to 
126 cycles of the frequency to = 1.7 eV and T w 300 fs. 

Although the details of the model potential do not play a significant role, at least on a 
qualitative scale, it is useful to know the level scheme in order to understand resonances 
observed in the numerically obtained harmonic spectra. With the chosen parameters the 
lowest four excited states have energies Qi = 0.985, Q 2 = 1-375, ^3 = 1.548, and ^4 = 1.592. 

In Fig.|2] we present examples of harmonic spectra obtained by Fourier-transforming in 
time the expectation value (x(t)) = J J dpdip p^*(p,ip, t)p cos </? \& (pip, t) and plotting the 
square of the result (hereafter called dipole "strength"). As expected, the structure of the 
spectra follows the selection rule, confirming the filtering effect, i.e., the number of harmonics 
present in a certain frequency range decreases with increasing rj. Additional lines of small 
intensity are also present. Those lines are due either to decays from excited states which are 
populated during the laser pulse or to recombinations of the electron with states different 
from the ground state. Such phenomena are also present in the interaction between a laser 
field and a circular molecule and are extensively discussed in [0. These extra lines can be 
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particularly useful for deriving informations about energy shifts because of the dynamical 
Stark effect. 

In order to verify the analytical findings discussed in the previous section a series of nu- 
merical simulations have been performed. In particular, given a fixed value for E 2 , a series of 
simulations have been run for different E\ and 77. In Fig.|3| the highest resolvable harmonic 
obtained from the simulations is compared with the highest harmonic expected from the 
calculations. The interpretation of Fig.||| is not straightforward and requires some discus- 
sion. Because of angular momentum conservation, in the two-color scheme with opposite 
polarizations the emission of a harmonic is possible only when \K W — K V uj\ = 1, where K w 
and K V u are the number of photons absorbed from the first and the second laser, respec- 
tively. Therefore, in order to achieve an efficient harmonic emission, it is required to find 
a regime of frequencies and field intensities where the absorption of K u and K vw photons, 
respectively, has a reasonably high probability. It follows that the extension of the harmonic 
spectrum, as predicted by the analytical calculation, plays the role of an upper limit. In fact, 
the analytical calculation in the previous section does not incorporate the actual absorption 
processes and their amplitudes, but the electron is rather "put by hand" into the continuum. 
This is the same for the well known case of a single linearly polarized laser where we have 
the predicted cut-off at I p + 3.17U P . However, this is verified only if the laser frequency and 
intensity are chosen within proper ranges so that the laser frequency is much smaller than 
I p , and the ponderomotive energy is comparable or larger than I p . 

Finding the combination of laser parameters that yields the most efficient harmonic gen- 
eration is not straightforward. However, a first hint about the most promising region in 
parameter space can be obtained by considering the absorption processes from the two 
lasers as independent. The three plots of the cut-off vs. field amplitude E\, shown in Fig.|| 
were calculated for fixed E 2 = 0.16 and different frequency ratios 77 = 3, 4, 5. The agreement 
between numerical simulations and analytical calculations is particularly good for 77 = 5 
but less satisfactory for 77 = 3 and 77 = 4 where the extension of the numerically obtained 
spectrum is less than the expected cut-off. By choosing other laser intensities it is possible 
to have the good agreement for, e.g., rj — 3 instead of 77 = 5. However, the set of plots in 
Fig. [3] demonstrates that it is possible to obtain harmonic spectra with significant extension 
even though the laser parameters are not precisely optimized. Note, that in Fig. |3]c the 
highest observable harmonic in the numerically obtained spectra is slightly greater than the 
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predicted cut-off. This is also well-known from the linearly polarized case and can be at- 
tributed to the fact that in the analytical calculations it is assumed that the electron is born 
and recombines exactly at the origin (where the nucleus is located). Allowing for offsets 
from the origin also yields harmonics beyond the calculated cut-offs. 

Another feature that should be noted is the asymptotic behavior of the highest observable 
harmonic order in the numerical simulations. Increasing the electric field amplitude Ex, 
the probability of absorbing a certain number of photons K from the first laser increases. 
However, the emission of high harmonics requires not only the absorption of many photons 
K from that laser but also the absorption of K±l photons from the other laser. This second 
part of the process is the real constraint. In fact, in Fig. [3] the value of E 2 is always constant 
and the probability of absorbing K photons decreases very rapidly with K and cannot be 
compensated with the increase of the probability of absorbing K ± 1 photons from the first 
laser. The same effect can be observed by looking at low values of E\ in Fig. ^. Here it 
is the low probability of absorbing photons of frequency uo which suppresses the harmonic 
generation. Summarizing, one can state that the agreement between theory and simulations 
is good when the absorption processes of order (I p + / y ri U)/(r] + 1) ± 1 have a reasonably 
high probability for both laser fields. 

The intensity of the emitted harmonics plays, of course, a key role in harmonic generation. 
Finding a configuration which enhances the efficiency of harmonic generation is important 
for possible applications because it may allow to use less intense lasers for obtaining a 
desired radiation intensity. In Fig. f| the spectra obtained with three different values of r] 
are compared. The harmonics in the rj = 4-spectrum are significantly more intense than 
those obtained with rj = 2 over a wide range of harmonic orders, although the scaling of 
the theoretical cut-off, i.e., a decreasing cut-off order with increasing 77, may suggest that 
the opposite should be true. However, our version of the Lewenstein model yields only an 
upper limit for the cut-off. The fact that higher 77 is favorable here can be understood 
considering that given a certain harmonic order, the number of photons required for the 
emission is inverse proportional to (77 + 1), and a lower order absorption process is (for the 
laser parameters chosen) more likely than a higher one. 

So far, only harmonic spectra for a given value of the electric field E 2 have been discussed. 
In Fig. |5] we present three different spectra obtained for different fields E 2 and keeping Ex 
constant. When E 2 is low only a few harmonics are present. With increasing field E 2 the 
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spectrum assumes a plateau structure. Increasing E% further leads to violent ionization and, 
thus, inefficient harmonic generation. 

Another feature predicted by the analytical calculations in the previous section is the 
intensity of the two harmonics g(rj + 1) ± 1 of a couple g with respect to the two laser 
intensities. From Eq. and (^) follows that the intensities of the two harmonics should 
be very different for E\ <C E 2 and should become very close for E\ ^> E 2 . In Fig.|B| the 
dipole strengths of the two harmonics in the couple g = 2 are plotted versus E\ for 77 = 3 
and 7] = 5. The numerical result clearly confirms the expected behavior. All the harmonics 
in Fig. ^| are obtained through the absorption of two photons from the second laser. We have 
chosen those couples because the total absorbed energy is below the ionization energy and 
resonance or interference phenomena play a marginal role. 

At low intensities it is possible to consider the absorption from the two lasers as two 
distinct processes. Therefore we can write the probability T(K W , K VUJ ) for the overall process 
as 

r(K ut K vu ) = r(K u )r(K flu ). (27) 

In terms of the generalized cross sections gk w and ctr: one has 

Y{K m K^) oc a Ku El K «a Knw El K «\ (28) 

Note, that in contrast to the single, linearly polarized laser field where the nth harmonic 
perturbatively scales with the nth power of the laser intensity /, here, in the two-color case, 
this is not true. Given for example the [g(rj + 1) + l]th harmonic, it scales with the (g + l)th 
power of the intensity I± and the gth power of the intensity I 2 . Consequently, by plotting 
double-logarithmically the intensity of a certain harmonic versus E\ or E 2 one expects a 
straight line, the slope of which yields information about the number of photons absorbed 
from that laser. Instead, if no straight line is obtained, the factorization of Eq. (pTp is not 
valid. In Fig. [7] the low-intensity region of Fig.|6|b is plotted on a log- log scale. We see that all 



the points are indeed aligned along straight lines so that (j27| ) is an acceptable approximation 
there. However, moving towards higher field amplitudes E\ in Fig.|| it is obvious that the 
approximation ([27]) will soon break down. 

Finally, let us briefly discuss the case where the two laser fields have the same circular 
polarizations and the selection rule g(rj — 1) ± 1 holds 0. In order to achieve angular 
momentum conservation in the emission process of a certain harmonic, the absorption of 
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photons from one field has to be accompanied by the emission of photons of frequency equal 
to that one of the other laser field ||10| . According to the selection rule, a frequency ratio 
rj = 5 in the co-rotating configuration provides a spectrum with the same harmonic orders 
present as with rj = 3 for counter-rotating electric field vectors. The two spectra are shown 
in Fig. |^. We see that in the region around harmonic order n = 50 (off all resonances) the 
intensity of the harmonics in the co-rotating case is significantly lower (between one and two 
orders of magnitude). Given a certain harmonic, the relative intensities can provide useful 
information for further investigations about the absorption and emission processes. 

IV. CONCLUSION 

We have investigated the generation of harmonics by atoms or ions in the two-color, 
coplanar field configuration for different values of integer frequency ratio and different laser 
intensities. Through an analytical calculation based on the Lewenstein model, the selection 
rule for the harmonic orders in this field configuration, a generalized cut-off for the harmonic 
spectra, and an integral expression for the harmonic dipole strength has been calculated. 

Numerical ab initio simulations of a two-dimensional model ion subject to the two-color, 
coplanar field configuration were performed. The numerical results did not suffer from the 
various assumptions made in the Lewenstein model and therefore served as an important 
benchmark for the theoretical predictions. 

The scaling of the cut-off as a function of both, one of the laser intensities and frequency 
ratio rj, as well as entire spectra for different t] and laser intensities were presented and 
analyzed. The theoretical cut-off was found to be an upper limit for the numerical results. 
The theoretically predicted relative strength of the two harmonics g(r] + 1) ± 1 in a certain 
couple g = 1, 2, 3, . . . when one laser is much more intense than the other was confirmed by 
the numerical simulations. The dipole strength of the harmonics in general increase with 
increasing r\ although the scaling of the theoretical cut-off, i.e., a decreasing cut-off order 
with increasing 77, may suggest the opposite. This was found to be due to the decreasing 
order of the absorption processes involved. 

Further studies can be undertaken in order to investigate the importance of resonances 
within this scheme and the possibility of using such resonances for enhancing strongly the 
efficiency of certain harmonic lines. 
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FIG. 1: The function C\{t)/U is plotted for rj = 1,2,3,4. The time r is counted in periods 2tt/oj. 
For rj = 1 the maximum energy gain is, as expected, 3.17 U p , for higher rj the maximum energy 
decreases. 
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FIG. 2: Harmonic spectra for r\ = 3 and r\ = 5. Both spectra have the structure expected from the 
selection rules. Additional lines between the 15th and the 30th harmonic are due to the population 
of excited states. The electric fields are E\ = 0.16 and E<i = 0.13. 
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FIG. 3: Highest observable harmonic as expected from the Lewenstein-type model (*) and from 
the numerical simulations (+) for different values of r/ versus the electric field E\. The electric 
field E2 is constant, E2 = 0.16. With this particular laser intensity the agreement between the two 
curves is good for rj = 5. 
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FIG. 4: Harmonic spectra for different values of rj. Diamonds: r\ = 2, triangles: r/ = 3, squares: 
i] = 4. The laser fields were E\ = 0.13 and E<i = 0.16. In all three cases the harmonic spectra 
show a similar structure. The conversion efficiency increases with rj. 
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FIG. 5: Harmonic spectra for different values of the electric field E^. Diamonds: E% = 0.06, 
triangles: E2 = 0.13, squares: E% = 0.18. Th electric field E\ is constant, E\ = 0.13. With 
increasing value of the electric field Ei a "plateau-like" structure in the spectrum appears. 
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FIG. 6: Behavior of the two harmonics belonging to the second couple vs. the laser intensity E\. 
E2 = 0.13 was held constant. While the 7th harmonic for rj = 3 and the 11th harmonic for r/ = 5 
are given by the absorption of two photons from the second laser and one from the first laser, 
the 9th harmonic for r/ = 3 and the 13th harmonic for r\ = 5 are generated by the absorption of 
two photons from the second laser and three from the first one. As expected from the Lewenstein 
model, with increasing intensity of the electric field E\, the strengths of the two lines of each couple 
become closer. 
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FIG. 7: Dipole strength of the two harmonics no. 9 and 11 of the second couple (rj = 4) vs. the 
electric field E\. The amplitude E2 = 0.13 was held constant. All the points are aligned along a 
straight line on the log-log scale. The slopes are 2 and 6, corresponding to the absorption of one 
and three photons, respectively. 



20 



10° 




• 



-25 



20 40 

Harmonic order 



60 



80 



FIG. 8: Harmonic spectra for opposite and same polarization of the two laser fields with rj adjusted 
in such a way that the same selection rule 4g ± 1 holds. With rj = 3 (diamonds) the two laser 
fields have opposite polarization, while for r\ = 5 (triangles) the polarization is the same. While 
the emitted harmonics are the same, the intensity (apart effects due to resonances) is significantly 
higher in the case where the laser fields have opposite polarization. The electric fields are E\ = 0.16 
and E 2 = 0.13. 
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